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Abstract. The rapid experimental progress in the field of superconducting nanocir- 
cuits gives rise to an increasing quest for advanced quantum-control techniques 
for these macroscopically coherent systems. Here we demonstrate theoretically 
that stimulated Raman adiabatic passage (STIRAP), a well-established method 
in quantum optics, should be possible with the quantronium setup of a Cooper- 
pair box. We find the parameters which optimize the procedure and show how the 
scheme appears to be robust against decoherence and should be realizable even 
with the existing technology. 



1 Introduction 

During the past decade the successful development of nanotechnology has made it possible to 
observe macroscopic quantum phenomena in superconducting nanocircuits based on Joseph- 
son junctions. These circuits behave like "artificial atoms" , and phenomena, traditionally part 
of nuclear magnetic resonance (NMR), quantum optics and cavity quantum electrodynamics, 
have been successfully implemented in such systems. Rabi oscillations and Ramsey interferom- 
etry [1|2|3| . the detection of geometric phases [4], two coupled artificial atoms |5|6j . artificial 
atoms coupled to electromagnetic resonators [718] , cooling techniques [9] , an analogue of elec- 
tromagnetically induced transparency [lOj and adiabatic passage in superconducting nanocir- 
cuits |ll|12|13|14|15j are examples first realized in atomic systems that have also been recently 
demonstrated with superconducting quantum circuits. 

One of the strongest motivations of these studies is to apply more complex control techniques 
of the dynamics in a nanodevice. Advanced control has been proposed as a tool to achieve 
dynamical decoupling from solid-state noise [16] . Compared to quantum optics and NMR, in 
these systems coupling energies larger than in atomic physics are easily achieved, thus reducing 
typical time scales for operations. Flexibility in the design offers several solutions for tuning 
couplings, allowing in principle to implement the Hamiltonian and to achieve any desired state 
transformation. On the other hand, the macroscopic nature of the system and the effects of 
solid state noise impose obstacles requiring more than a mere translation from quantum optics 
to solid state. 

A particularly interesting technique in quantum optics is the so-called stimulated Raman 
adiabatic passage (STIRAP), developed by Bergmann and co-workers [17]. This method rep- 
resents one of the most efficient coherent population transfer schemes between two (or more) 
quantum states known in quantum optics. It has many applications in such diverse areas as 
chemical-reaction |18j . laser-induced cooling [19] , single photon generation in atom-cavity sys- 
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tern |20|21j and represents a building block of non-abelian geometric quantum computation 
procedures [22]. 

In this work we will demonstrate the application of the STIRAP technique to a single 
Cooper-pair box in the charge-phase regime (the so-called quantronium) [2]. Although one 
can apply this procedure to different regimes and setups of superconducting nanocircuits 
[14I15I12I13) . we have chosen the quantronium because it is very similar to the atom-laser 
system in quantum optics and it has been extensively studied with respect to its decoher- 
ence properties. Actually the main problem in nanodevices is the trade-off between efficient 
coupling to the driving field and protection against noise. In this respect quantronium offers 
a uniquely convenient design, where tunability is exploited to obtain selective and relatively 
strong coupling to the fields allowing to perform STIRAP before decoherence takes place. 

In Section 2 we will describe the STIRAP protocol in the quantronium. In Section 3 we 
will find the parameters which optimize the efficiency of the population transfer and finally in 
Section 4 we will discuss the feasibility of the protocol in the presence of decoherence. 



2 STIRAP in the quantronium 

The STIRAP technique in quantum optics is based on a three state linkage forming a A pattern 
of two hyperfine ground states \g) and \u) coupled to an excited state |e) by classical laser fields 
A g (t) coscugt, A u (t) cosuj u t |17|23j . Applying the unitary transformation diag(l, e 1 "" 4 , e luJgt ) in 
a basis {|e), \u), \g}} and the rotating wave approximation, the three-level Hamiltonian reads 
(ft = 1) 

H = 5 s \e)(e\ + (6 s -5 u )\u)(u\ + ±(A u (t)\e)(u\+A s (t)\e)(g\+h.c.) (1) 

where S s = E c — E g — o; g and <5 U — E e — E u — lu u are the detunings and A u (t), A s (t) the 
slowly-varying Rabi frequencies. 

An essential condition for STIRAP is that there be two photon-resonance between states 
\g) and \u) meaning 

5 = S g - 6 U = (2) 
Under this condition, setting A = 8 g = S u , the Hamiltonian becomes 

H = A\e)(e\++±(A u (t)\e)(u\+A e (t)\e)(g\+h.c.) (3) 

The peculiarity of this Hamiltonian is the presence of a null instantaneous eigenvalue, whose 
associated adiabatic state is essential for the population transfer mechanism. It has no compo- 
nent of the excited state |e), being a coherent superposition of the states \g) and \u) only and 
can be written as 

\<j>°) = \g) cos6>0) - \u) sin0(t) (4) 
where the time-dependent mixing angle is defined as 

*».(,) = Mj (5) 

The latter state cannot decay by spontaneous emission from e) and is usually called dark state. 
From Eq. ((4]) it can be seen that by slowly varying the coupling amplitudes A u (t), A s (t) the 
dark state can be rotated in the two dimensional subspace spanned by \g) and \u). STIRAP is 
based on tying the state vector to the dark state. Population transfer from the state \g) to the 
state \u) can be achieved by applying the so-called counterintuitive scheme |17j : the system is 
prepared in the state \g) at vanishing couplings and the state vector coincides with the dark 
state; then the field ^4 U is turned on while A g is kept zero. Finally, by slowly switching ^4 U 
off while A g is switched on, the dark state evolves from the bare state \g) to the target state 
\u). If the adiabatic condition is fulfilled, the state vector will "follow" the dark state and the 
population transfer is achieved. 
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Fig. 1. Idealized circuit diagram of the quantronium. This quantum-coherent circuit consists of a 
Cooper pair box island (black node) delimited by two small Josephson junctions (crossed boxes) in a 
superconducting loop. The bigger junction with Josephson energy £j is used for readout. The offset 
charge of the island can be tuned with the gate voltage Vg (C g >C C), while the Josephson energy of 
the junctions can be controlled via the magnetic flux <£>. Here we choose 5& = and Ej/Ec = 1.318. 



Our aim now is to implement the Hamiltonian ([3]) in the quantronium and realize adiabatic 
population transfer. The quantronium circuit [2] is showed in Fig. [1] It consists in a supercon- 
ducting loop interrupted by two adjacent tunnel junctions with Josephson energies Ej/2 and by 
the larger readout junction (with ej Ej). The two small junctions define the superconducting 
island of the box, whose total capacitance is C and charging energy Eq = (2e) 2 /2C. Except at 
readout, the circuit can be modelled in the laboratory frame by the following Hamiltonian 

H = Y J Ec[n- n g {t)] 2 \n)(n\ - ^-(\n)(n + 1| + h.c.) (6) 

n 

where {|n)} are eigenstates of the number operator N of extra Cooper pairs on the island and 
rig(i) = C g V g /(2e) is the reduced gate charge. The Cooper pair number N and the supercon- 
ducting phase difference across the readout junction are the degrees of freedom of the circuit. 
A dc voltage V g applied to the gate capacitance C g and a dc current applied to a coil producing 
a flux in the circuit loop, tune the quantum energy levels. Here we choose = and a ratio 
Ej/E c = 1.3 {E c = 0MK B K, Ej = 0.87K B K as reported in Ref. [2]). Being Ej ~ E c , neither 
N nor the phase different is a good quantum number. The box thus has discrete quantum states 
which are quantum superpositions of several charge states with different n. Diagonalizing the 
Hamiltonian ((6]) for a certain dc bias voltage n g (t) — n g o, one finds 

i 

where {|</>i)} are the eigenstates of the undriven system. In Fig. [5] are plotted the lowest four 
energy levels of the system as a function of the reduced gate charge n g . At the particular bias 
point n g o = 1/2, the system is immune to first-order fluctuations of the gate charge and the 
lowest decoherence rates are obtained |2|27j . 

Manipulation of the quantum state is performed by adding to the dc part of the gate voltage 
ac microwave pulses with small amplitudes n g (t) — n g o+rig C (t). Adding the ac part of the drive, 
the Hamiltonian in the basis of the undriven system eigenstates, becomes 

H = Ei I <tH ) I + <(*) E n ^ I & ) I ( 8 ) 

where riij = (<f>i\N\<pj) are the matrix elements of the number operator in the basis of the 
undriven eigenstates. 

The STIRAP protocol can be carried out between the three lowest energy levels with an ac 
gate charge n a g c = A g (t) cos uj g t + A u (t) cosui u t (see FigJ5]). Besides the minimal decoherence, 
at the optimal working point n g o = 1/2, the microwave field couples off-diagonally in the basis 
of the eigenstates. However at this point small level spacings and selection rules impede the 
operation of the scheme [15] . This problem can be circumvented by working slightly away from 
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Fig. 2. The four lowest energy levels of the quantronium as a function of the gate charge n g . At the 
chosen working point n g o, STIRAP can be performed between the three lowest energy levels \g), \e) 
and \u) with resonance frequencies Tiu) g = E e — E g and hui u — E e — E u . 



the optimal point, but not too far to avoid high decoherence. We have chosen n g o = 0.47 (see 
Section 13.11 for more details) . Here the Hamiltonian is more complicated than the ideal one of 
Eq. ([3|), because of the presence of more than three levels, diagonal and counter-rotating drives. 
Furthermore, the effective peak Rabi frequencies are different (see Sec. I3.3[) and the coupling 
between the states \g) and \u) is not zero (see Fig. 0k). Nevertheless, if the two ac signals with 
slightly detuned frequencies u> g and u> u are applied to the gate, it is possible to adiabatically 
transfer the population from the ground state \g) to the first excited state \u) just as in a 
quantum optical system. 

The von Neumann equation for the isolated system is 

P = \\p,S\ (9) 

where p is the density matrix of the system. In Fig. [3j there we show the numerical solution 
of this equation for the Hamiltonian ([8]) truncated to eight charge states. First the system is 
prepared in the state \g). Then, two Gaussian-shaped microwave pulses are applied in the coun- 
terintuitive sequence and finally a population transfer with almost unit efficiency is obtained. 
The state |e) remains basically unpopulated during the whole procedure as one would expect 
from the STIRAP protocol. 




-7500 -5000 -2500 2500 5000 7500 -7500 -5000 -2500 2500 5000 7500 

time t time t 



Fig. 3. Population transfer by STIRAP in the quantronium. a) Gaussian-shaped microwave pulses 
A s (t) = Aexp [~(t - t) 2 /T 2 ], A u (t) = Aexp [-(t + r) 2 /T 2 ] with a maximum amplitude A = 0.05. For 
a charging energy of Ec = 56.9/ieV the time unit corresponds to about 1.16 x 10 _11 s. The pulse width 
is T — 20ns and the delay between the two pulses r = 14ns. Detunings of the frequencies u s and uj u 
are zero, b) Time evolution of the population p gg , p uu and p ee for the isolated system. The efficiency 
achieved for population transfer is 99.9%. 
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3 Sensitivity to parameters 



Once we have shown how the STIRAP procedure can be implemented in the quantronium, in 
this section we discuss the parameters which optimize the transfer efficiency. 



3.1 Optimal working point 



As described in the previous section, it would be preferable to work at the optimal point 
n g0 — 1/2, where decoherence is minimal and the drive is purely off-diagonal. However at 
this protected point, unwanted transitions to higher levels can occur and, moreover, the parity 
of the wavefuntions induces selection rules |15j . which render the implementation of STIRAP 
protocol impossible. In particular the selection rules make the coupling between the states \g) 
and |e) vanish at the optimal point (see Fig. 2^). The solution, then, is to work away from 
this point in order to have enough coupling between the above-mentioned states. On the other 
hand, moving excessively away from the optimal point leads to high decoherence rates and 
quasi-particle tunneling processes [27] , 

In Fig.0J), we plot the populations p uu , p cc and p gg at the end of the STIRAP procedure as a 
function of the bias gate charge n g Q. One can immediately see that at n g o = 1/2 the population 
transfer is not achieved and the system remains in the initial state \g). Moreover, the further 
we move away from the protected point, the more the transfer efficiency (which coincides with 
the population p uu ) grows due to the increase of the coupling n gc . A good choice which assures 
high efficiency and acceptable decoherence rates is n g o = 0.47. 
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Fig. 4. a) Charge matrix elements for the three levels involved in the STIRAP operation. At the point 
n g o = 1/2, the coupling term n sc — due to selection rules, b) Populations p uu , p Ge and p gg at the end 
of the procedure as a function of the gate charge. A pulse width T = 20ns, a pulse delay r = 14ns and 
Gaussian microwave pulses with maximum amplitude A = 0.05 have been used. In absence of noise 
and with the parameters chosen, population transfer is efficient for n g o < 0.47. 



3.2 Sensitivity to Delay 



The following considerations concern the choice of an optimal delay between the two pulses. 
The optimal delay, which leads to maximal transfer efficiency, must be chosen by maximizing 
adiabaticity. The adiabatic condition requires \ Aj/Aj\ < ujj (j — u,g). This condition translates 
in an upper bound on the time derivative of the mixing angle 0(t), which is proportional to 
non-adiabatic couplings [17] . Thus for maximal adiabaticity the mixing angle 9(t) must change 
slowly and smoothly in time in order to guarantee small non-adiabatic couplings . 

The numerical simulation in Fig.[5]shows that for our choice of Gaussian pulses, the optimum 
occurs when the delay is slightly lower than the pulse width. 
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Fig. 5. Transfer efficiency as a function of the pulse delay. The optimal delay is slightly lower than the 
pulse width T ~ 20ns for Gaussian-shaped pulses. In the low efficiency regions non-adiabatic couplings 
between the instantaneous eigenstates occur. 

3.3 Sensitivity to Rabi frequencies 

For STIRAP it is better to have two nearly equal peak Rabi frequencies (maximum drive 
amplitudes). But looking at Fig. [4^l, we see that the charge matrix element coupling the states 
\g) and |e) (n ge ) is smaller than the one coupling the states |e) and \u) (n eu ). This means 
that the effective peak Rabi frequencies are different n ge A g < n U0 A u . It is interesting to note 
that, in our analogy with quantum optics, the charge matrix elements correspond to the dipole 
transition moments. 

If the two maximum Rabi frequencies are different, while the pulse widths are about the 
same, the projection of the state vector onto the adiabatic transfer state is very good initially 
(because in our case the more intense drive comes first), but necessarily less good in the final 
stage and consequently the transfer efficiency will be small. 

It is also true that using large pulse areas, small deviations from the optimal conditions do 
not lead to significant drop in transfer efficiency. In fact adiabaticity requires a large pulse area 
[17j and non-adiabatic transitions can be suppressed by increasing the drive amplitudes. How- 
ever using equal peak Rabi frequencies and equal pulse widths allows to reduce the necessary 
pulse areas and facilitate efficient population transfer. Furthermore in solid state qubits it is not 
easy to increase the microwave pulse amplitude beyond a certain value since other excitations 
come into play. Thus it is important to externally adjust the drive amplitudes in such a way 
that the effective Rabi frequencies are about the same (n ge ^4 g ~ n ue A u ). 

In Fig. [6j we report a numerical plot of the transfer efficiency as a function of the ratio of 
the two peak Rabi frequencies. Simply making the drive A g twice as strong as the drive A Ul we 
get a high transfer efficiency with reasonable amplitudes. 

3.4 Sensitivity to detunings 

When the two frequencies u g and ui u are not exactly resonant with the respective transitions, 
the presence of non-zero detunings can affect the efficiency. There are two cases of interest. First, 
both uj g and w u may vary, while the two-photon resonant condition is maintained (single-photon 
detuning, 5 = 0). Alternatively either frequency may vary, while the other is fixed (two-photon 
detuning, 5^0). 

The single-photon detuning does not affect the formation of the dark state, because the 
mixing angle does not depend on it. However this type of detuning affects the adiabatic con- 
dition [17] and when we increase i5 g (and consequently 5 U ), the transfer efficiency decreases for 
deterioration of adiabaticity. 

The detuning from two-photon resonance is more detrimental for STIRAP, because it pre- 
vents the exclusive population of the trapped state, which is no longer an instantaneous eigen- 
state of the Hamiltonian. A more detailed analysis of the instantaneous eigenstates when 5^0 
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Fig. 6. Transfer efficiency as a function of the ratio of the maximum drive amplitudes. The simulation 
has been done with a pulse width T — 20ns, a pulse delay r = 14ns. The ratio of the Rabi frequencies 
has been changed by keeping fixed the amplitude A u — 0.02 and varying A s . The value A e /A u — 2 
assures a high transfer efficiency . 



shows that there is no adiabatic transfer state providing an adiabatic connection from the initial 
to the target state, as does the dark state for S = and the evolution leads to complete popu- 
lation return of the system to its initial state. The only mechanism which leads to population 
transfer is by non-adiabatic transitions between the adiabatic states [17]. Actually for small 
values of 6, narrow avoided crossings between the instantaneous eigenvalues can occur and the 
population can be transferred by Landau-Zener tunneling [24) . 

To show how STIRAP efficiency is much less sensitive to single photon detuning than to two- 
photon detuning, we show in Fig. [7]the numerically calculated transfer efficiency plotted versus 
the single-photon and the two-photon detuning. The figure shows two plots of the efficiency 
sensitivity to detunings for equal drive amplitudes (Fig. [7^,) and for externally adjusted drive 
amplitudes (Fig [7b). In the case of adjusted peak Rabi frequencies (see section l373|) . the region 
of high efficiency is larger. This is important for reducing the effect of low frequency noise 
treated in section |4~21 



1 




8 - 8 U (units of T-') 5 - 5 u (unitsofT" 1 ) 



Fig. 7. Transfer efficiency as a function of single-photon and two-photon detuning for equal peak Rabi 
frequencies a) (A s — A u — 0.03) and for externally adjusted amplitudes b) (A s = 2A U = 0.06). The 
region of high efficiency is larger in the figure b) showing the importance of choosing appropriately 
adjusted external amplitudes. 
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4 Effects of decoherence 

The efficiency of solid-state quantum-coherent nanodevices is affected by various (device depen- 
dent) noise sources. In the quantronium, high-frequency noise and low-frequency noise coexist. 
The former is mainly responsible for unwanted transitions which lead to exponential decay of 
the signal in time, whereas the latter mainly makes the calibration of the device unstable and 
determines the initial power-law decay of the signal |25|27j . 

Based on the above considerations, we discuss the feasibility of the protocol without making 
a detailed analysis of decoherence in the protocol [26] , which is beyond the scope of this work. 
The key observation is that strong unwanted processes involving the state |e) have negligible 
effects both for high and for low-frequency noise on the STIRAP protocol. Decoherence is mainly 
determined by processes involving the states \g) and which have been well characterized in 
the quantronium and, as a matter of fact, allow decoherence times tr > 300 nsp7j. 

4.1 High frequency noise 

High frequency noise can be analyzed by the quantum-optical master equation 

P = ~[ P ,H']-r P (io) 

where p is the density matrix of the system and H' is the Hamiltonian ([5J) in the rotating frame 
[28) . The dissipator Ep includes spontaneous decay rates and environment-assisted absorption 
in the presence of the driving fields. In the basis {\g), |e), u}}, it reads 

( r P)ij = 1% 2 ^ Pi 3 ~ (! - S ij)lPij ~ % Pkklk^i (11) 

k 

where 7$ = X)fc=£j Yi— >k- We take the dissipator time-independent (which overestimates decoher- 
ence) and include all transitions as well as dephasing rate 7, which accounts phenomenologically 
for low- frequency noise. For the state |e) (the second excited state of the quantronium) we as- 
sume 7c = 7o^ u + 7c^g = 27 u where 7 U is taken on the order of the rate observed in the 
experiments of Ref. [27]. While in quantum optics the rate 7 u ^ g vanishes and the remaining 
decay rates act on depopulated states hardly affecting the protocol, STIRAP in the quantro- 
nium may be sensitive to the extra decay \u) — ► \g) involving the two low-lying states. The 
Fig. [8] show results for the master equation (|10| . It can be recognized a good robustness of the 
procedure against decoherence. We observe that the extra rate 7 u ^ g changes the population 
only during the waiting time after completion of the pulse sequence and the population of level 
e) is slightly increased. 

4.2 Low frequency noise: detuning instability 

Low frequency noise is modelled as due to impurities which can be considered static during each 
run of the protocol but switch on a longer time scale, thus leading to a statistic distribution 
of level separations. In protocols like Ramsey interference, they result in a distribution of 
oscillating frequencies, whose average determines defocusing of the signal [25]. The effect of 
averaging over static impurities configurations corresponds to statistically distributed energy 
level fluctuations 5E\, which, in the rotating frame, translates in fluctuations of the detunings. 
As one can see from Fig. [9] even large fluctuations of E e (X2) hardly affect STIRAP, since they 
still leave equal detunings of both microwave fields (i.e. they do not affect two-photon resonance, 
see section l3"/f]) . On the other hand, fluctuations of the difference between the energies of the 
two lowest eigenstates E n — E g (XI) are potentially detrimental since they translate in two- 
photon detuning fluctuations (see section l3~4]) . However as long as the system is well-protected 
from noise, fluctuations remain inside the high efficiency region in Fig. [7] and almost complete 
population transfer can be achieved. 
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Fig. 8. Time evolution of the populations p gg , p ee and p uu with decoherence. The time unit corresponds 
to about 1.16 x 10~ n s {h/Ec). For the calculations with quantum noise we have used the decay rate 
7 U = 4.4 x 10~ 5 and the dephasing rate 7 = 2.6 x 10 -4 |27] . The other parameters are those used in 
Fig. [3] The extra decay 7 u ^ g acts only at the end of the procedure and the increase in the population 
of level |e) is harmless. 
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0.2 0.47 0.8 1 

Fig. 9. Detuning instability due to low-frequency noise. Fluctuations in the gate charge An g translate 
in fluctuating single-photon detuning S g — > S s + X2 and fluctuating two-photon detuning 8-^5 + Xi . 

5 Conclusions 

In summary we have discussed the possibility to implement STIRAP in the quantronium device. 
One important advantage of the protocol is that the efficiency does not depend sensitively on 
details of the procedure and on timing making it robust against moderate fluctuations of the 
solid-state environment. One goal of this work is to provide a simple basis and parameters for 
the experimental realization of STIRAP. Since we have used real circuit design with common 
parameters and corresponding time scale, it should be feasible to experimentally verify STIRAP 
in nanocircuits with the state-of-the-art technology. 

There are many interesting applications of this scheme such as the preparation of Fock 
states and the single photon generation [28] in cavity coupled to the nanocircuit. The cavity 
can be implemented using electrical resonator [29], transmission lines [7] and nanomechanical 
resonators [T3] . 
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